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Abstract 



Top-down models for the origin of ultra high energy cosmic rays (UHECR's) 
propose that these events are the decay products of relic superheavy metastable 
particles, usually called X particles. These particles can be produced in the re- 
heating period following the inflationary epoch of the early Universe. We obtain 
constraints on some parameters such as the lifetime and direct couplings of the 
X-particle to the inflaton field from the requirement that they are responsible 
for the observed UHECR flux. 
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1 Introduction 



The ultra high energy cosmic rays (UHECR's), events with primary energies 
above E > 10^" eV, constitute one of the most intriguing mysteries in nature. 
Their origin, source and composition are not understood since they were first 
observed by the Volcano Ranch experiment in f963 0. 

According to conventional theories, cosmic rays are produced by the accel- 
eration of charged particles in intense electromagnetic fields of astrophysical 
objects. This conventional wisdom, however, has some difficulties in explaining 
the observed flux of UHECR's. Most of the known astrophysical objects do 
not have an electromagnetic field strong enough neither a large acceleration site 
capable to accelerate particles up to fO^^ eV 

Another problem is the propagation of the UHECR's through the intergalac- 
tic medium. If the particles are protons or nuclei, they should quickly degrade 
their energy due to interactions with the cosmic microwave background (CMB), 
producing a cut in the spectrum at energies above 3 x 10^ eV, the so-called GZK 
limit §-0. Simulations Q show that in order that protons with energy larger 
than 10^° eV arrive at Earth they must have traveled less than wlOOMpc, the 
attenuation length. A cut off is also expected if the primary particles are high 
energy photons in the range of 3 x fO^^ eV to fO^^ eV, which have attenuation 
length smaller than fOO Mpc, due to their interaction with various background 
photons (infrared/optical, CMB or universal radio background), mainly by pair 
production and inverse Compton scattering. These consecutive interactions pro- 
duce electromagnetic cascades with decreasing energy that end up below fOO 
GeV resulting in a gamma-ray background. 

Furthermore, UHECR's are usually not deflected by intergalactic magnetic 
fields, so their incoming directions should point to the sources. However, taking 
into account the arrival directions of the UHECR events ||] and constraining 
their origin to be in a distance of fOO Mpc around the Earth, no candidate 
sources were found so far (an exception is a suggestion that these events can be 
generated in M87 galaxy in the Virgo Cluster and made isotropic by the galactic 
magnetic field 0). 

A new scenario, the so-called top-down mechanism, was then proposed to 
explain these events [0. In this scenario, UHECR's would originate from the 
decay of some superheavy particles (mx 3> f 0^" eV) inside the GZK volume. 
They could be either emitted by topological defects at the present epoch 
or created in the pos-inflationary era jTsj. In the second case, where the super- 
massive particles must be metastable {tx > fO^" years), they could have been 
created perturbatively, by a mechanism that is called reheating [ p^ , or by no- 
perturbative process like preheating jTs] , instant preheating [ p^ or gravitational 
production . In this letter we will study bounds on the direct coupling of X- 
particles produced during the reheating phase to the inflaton field by requiring 
that they are sources of UHECR's. A preliminary version of this work appeared 
in 
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2 Observed UHECR's flux and X-particle abun- 
dance 



We will now derive a relation among the abundance fix , the lifetime tx and the 
mass mx of the X-particle in order to provide the observed flux of UHECR's 
on Earth. 

It is reasonable to assume that X-particles decay into quarks and leptons. 
Particles are produced following the quark-hadron fragmentation process de- 
scribed by QCD that results in the production of nucleons (« 10%) and pions 
(« 90%) distributed equally among the three charge states. Charged pions 
produce basically electrons and neutrinos and neutral pions decay into pho- 
tons. Photons dominate the primary spectrum over protons by a factor of w 6 
[ p^ -pof. Following Sigl and Bhattacharjee we will consider that the 

UHECR's are photons with energy E ~ 10^° eV and an attenuation length of 
1{E^) w lOMpc. It was recently suggested that photons can indeed provide 
the bulk of UHECR's without violating EGRET bounds on GeV gamma-ray 
backgrounds pl| |. 

The photon flux generated by the decay of an uniform distribution of X- 
particles is given by: 

where dN^/dE^ is the photon injection spectrum from X-decay and nx{= 
it-x/tx) is the X-decay rate. 

From the observed UHECR's flux |l|l, J{E ~ 10^ eV) = lO^^r (^^2 g Qgy^- 

^^/dE^ « E-- 



and using dN^/dE^ « E from a QCD model [O, one gets: 



n,=3.7xl0--f-^)'''%m-s-. (2) 



,1012 GeV 

So, we can obtain a relation between flxh^ and tx- 

rx=^-^^W(f^x/i') years, 3) 

nxmx \ mx J 

where V,x = px/pc, px is the X-particle mass density, pc — I0~^h'^ GeVcm"^ 
is the critical density, and h is the present value of Hubble's constant in units 
of 100kms~iMpc-i. 

It could be interesting that these unknown X-particles constitute a fraction 
of the cold dark matter (CDM). In this case, they could be concentrated in 
galactic haloes and should not be uniformly distributed along the whole universe, 
as assumed in (^. Considering just the contribution of our galactic halo the 
correct equation to describe the flux of UHECR's is: 
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where i?'* is the radius of the galactic halo and h\ is the decay rate of the 
X-particles clustered in the halo. 

The flux above can be written in terms of J-y, given in (|l|), so that = f J-y 
and: 



7" jyh 



= l5y.W^( ]( WlOMpcW p^^,, \ 

K^CDAihy VlOOKpcy V ^(^^7) / VO.SGeVcm-^; 

where Pqjjj^ is the cold dark matter energy density in the halo. We recall that 
if X-particles are a fraction of CDM, it can be written as p\ = sPcoM- "^^^ 
fraction e is the same everywhere in the universe, so that e = fix /i^cDM- 
The halo concentration of X-particles implies a modification of equation : 

Tx = 10'' fnxh' / lO^^GeV y/^ 
\ mx J 

Considering mx ~ 10^^ GeV, for a minimum tx of 10^'' years and a maximum 
^xh' of the order of 1, we obtain the following limits (for / w 10^) (see also 

My- 



for (f7x/i^) - 1 Tx 10^5 years, 

for Tx 10^° years ^ (f^x/i^) - 10 



(6) 



Therefore, cosmologically interesting abundances of X-particles which at the 
same time provide a solution to the UHECR events require a very large life-time 
Tx- These requirements are generally difficult to obtain in realistic theories, but 
examples do exist p2[. 



3 Production of X-particles in the reheating phase 

In the inflationary scenario, following the exponential expansion of the universe, 
particles are created by the oscillation of inflaton field about its minimum. 
The coupling of the inflaton with other fields becomes important and allows 
the transference of its energy to other particles. In the perturbative reheating 
process [Q, the infiaton decays mostly into relativistic particles (radiation) 
that reheat the universe. In this context, the superheavy particle X could 
be created by either radiation annihilation process (indirect production) or by 
direct inflaton decay — > XX (direct production), if such a coupling is allowed. 
The study of the latter mechanism is the main focus of this work. 
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In order to study quantitatively X-particle production during reheating, 
we have to solve the coupled differential Boltzmann equations for the energy 
densities of X particles (px), radiation (pr) and inflaton (p^) : 

+ 3iJp0 + = 0, (7) 

PR + AH PR - (1 - Bx)r^p^ - ^-^[{pxf - (pf )^] = 0, (8) 

nix 

PX + 3Hpx - "^BxT^p^ + ^^[{pxf - {p'^f] = 0, (9) 

mx 

where: 

H"^ = —{p^ + PR + px)- (10) 

In the above equations H is the expansion rate of the universe, Bx is the 
branching ratio for the — > XX decay, is the equilibrium value for the X- 
particle energy density and {cr\v\) is the thermal average of the X annihilation 
cross section times the M0ller flux factor. The third term in equation (ph has a 
factor mx/m^ that corrects a similar equation in [^. According to Q, the 
annihilation cross section can be estimated as {(t\v\) ~ o? jT"^ where g — \J A-Ka 
is the gauge coupling strength and T is the characteristic temperature for the 
process. V ^ is the inflaton total width and is obtained as a function of the 
reheating temperature, Trh-, assuming instantaneous conversion of the inflaton 
energy density into radiation: 

= 1.93 X 10-i« (^1^^ GeV. (11) 

It is important to notice that the reheating temperature (Trh) is not the highest 
temperature achieved by the universe during reheating. For example, in this 
paper we are using a simple chaotic inflation model where V — m?^(jp' /2 and 

= 10^3 Qgv. Following H it can be shown that, even if Trh = 10^ GeV, 
the universe can achieve temperatures as high as 10^^ GeV. In this case, for 
g = 1 we have {(t\v\) — lO^'^^ GeV^^, which allows the process 77 XX. 

It is convenient to work with co-moving and dimensionless variables and 
therefore we defined appropriate parameters and solved the equations nu- 
merically for the scaled variable X = pxmx~^a? , where a{t) is the scale factor. 
In this case, the new independent variable is a; = am^. 

In Figure 1, we illustrate the behavior of X{x). The curves A, B and C 
show two plateaus associated with two mechanisms of X-particle production: 
the first one is the indirect production (77 — > XX) and the second is the direct 
one ((/) XX). Modifying Bx and keeping Trh constant (B and C curves) 
only the second plateau changes, reflecting the dependence of the direct decay 
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on Bx- On the other hand, when Trh is much smaller than 10^ GeV, there is 
not enough energy for the indirect process and hence, the first plateau tends to 
disappear. 



4 Finding Qx 

We can find the X-particle mass density fix today using the present value 
{x = xq) of the radiation abundance {^nh^ = 4.3 x 10^^): 

^xh"^ ^ Px{xo) , , 

In thermal equilibrium, T cx g^^^^'^a^^ p^ , and hence pR can be written 

as: 

PrM = PRixRH) [^^-^ ) [^-^ ) (13) 

where is the effective number of degrees of freedom (we used g*{xRH) — 200). 

Assuming that X particles are created already decoupled, we can write that 
Pxixo) = Px{xrh){xrh/xoY, where: 

Px{xrh) = XFXRH~^m4,mx- (14) 

The parameter Xp is the X value in the second plateau (Figure 1) and 
corresponds to X{xrh)- So, we get: 

0^,2 o u2Px{xrh)Trh 

Pr{xrh) Jo 

n /'Trh\ Xprnxml 
= 1.5x10^ ^ (15) 

where nipi is the Planck mass and Tq — 2.37 x 10"^'^ GeV is the present-epoch 
CMB temperature. 



5 Numerical Results 

We worked with acceptable reheating temperatures in the range of lO^'^GeV < 
Trh > lO^GeV. The lower limit is imposed by nucleosyntesis and the upper 
one is chosen not to lead to overproduction of gravitinos ||25|| . The subscript i 
refers to the epoch before inflaton decay. For chaotic inflation models we have 
Hi « m^. As initial condition for the differential equations we choose Xi ~ 1. 

Solving the Boltzmann equations numerically for different Bx and Trh, we 
obtain the respective Xp and consequently its associated ilxh^, as it is shown 
in the table 1. 
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Bx 


Trh = 10'-' GeV 


Trh = 10^ GeV 


Trh = 10-' GeV 


Xf 




Xf 




Xf 


nxh^ 





3.3 X 10-" 


0.38 














1 X 10-1^ 


1.8 X 10-4 


20.97 


1.78 X 10-" 


2.06 X 10-5 


1.78 X 10--* 


2.06 X 10-^^ 


1 X 10-" 


1.78 


2.06 X 10^ 


1.78 


2.06 X 10-1 


1.78 


2.06 X 10-" 


1 X 10-^ 


1.78 X 10^ 


2.06 X 10' 


1.78 X 10^ 


2.06 X 10^ 


1.78 X 10^ 


2.06 X 10-* 


1 X 10-^ 


1.78 X 10" 


2.06 X 10" 


1.78 X 10" 


2.06 X 10^ 


1.78 X 10" 


2.06 


1 X 10-^ 


1.78 X 10« 


2.06 X 10^^ 


1.78 X 10** 


2.06 X 10^ 


1.78 X 10« 


2.06 X 10^ 



Table 1: Numerical results for the final abundance obtained by solving the 
Boltzmann equations, for Trh = 10^ 10^, 10-^ GeV and Bx < IQ-^. 



For Bx = 0, the X particle production is not effective if the reheating 
temperature is lower than T = 10^ GeV. From table 1 we also get the expressions 
below, related to the direct production process that prevails over the indirect 
one, for Bx > 10-^": 

Xf = 1.78 X 10"Sx, (16) 
Cixh^ =^ 2.06x10' (^^^Bx. (17) 

An excellent analytical approximation to the numerical results can be ob- 
tained by requiring that nx — n^Bx, in which case we find: 

XF = i-^Bx (18) 
87r mj, 



Substituting the numerical values the above equation reproduces (17) exactly. 

The maximum efficiency of the direct decay process {<f> —> XX) is obtained 
when Trh is minimum. Recently, low temperature reheating models were ex- 
plored in Taking for the minimum value of Trh — 10^^ GeV and using 
the upper limit for X abundance (ilxh^ — 1) we get the maximum value for 
the branching ratio of Bx — 10"^. 

More generally, we can find bounds on Bx by requiring that X-particles are 
sources of UHECR's, in which case the limits of (^ apply. These bounds are 
shown in figure 2 as a function of the reheating temperature for two limits on 
the abundance: flxh'^ < I {tx 10^^ years) and flxh'^ > 10"^^ {tx =i 10^° 
years). 



6 Conclusion 

We have investigated the constraints on the properties of supermassive relic 
metastable particles demanding that their decay products are responsible for 
the observed UHECR flux. In particular, we found the relevant range for the 
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branching ratio of inflaton decay into the X-particle as a function of the reheat- 
ing temperature. 

For very small branching ratio of the inflaton to the X-particle, Bx > 10~^^, 
the direct production (0 XX) dominates over the indirect thermal produc- 
tion even for Tr// = 10® GeV. For the smallest possible reheat temperature, 
Trh = 10~^ GeV, thermal production of X-particle does not occur and the 
direct mechanism with Bx < 10~^ is required. 

Further consequences of the top-down scenario are the dominance of photons 
and neutrinos as primaries of UHECR's (l9| and an anisotropy of these events 
as a result of the large concentration of X-particles in the galactic center [p7| . 
More recently, it was suggested that clumps of dark matter in the halo could 
explain the observed multiplets (12 doublets and 2 triplets) of UHECR's events 
[ p8| . A detailed analysis including several models for galactic haloes of our and 
other galaxies was performed in We look forward to new experimental 

results from current and future experiments to elucidate the composition and 
anisotropy of the UHECR events. In particular, the first site of the Auger 
Observatory in the southern hemisphere, being sensitive to the galactic center, 
will be important in testing the anisotropy resulting from the top-down scenario. 
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Figure 1: Evolution of X (= pxm^a^) as a function of the modified scale 
factor, X. We used to,/, = 10^^ GeV and mx = 10^^ GeV. 




Figure 2: Branching ratio Bx of the inflaton into X-particles as a function of 
the reheating temperature Trh for two Umiting values of flxh"^- 
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